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ABSTRACT 

We study the contamination of the B-mode of the Cosmic Microwave Background 
Polarization (CMBP) by Galactic synchrotron in the lowest emission regions of the 
sky. The 22.8-GHz polarization map of the 3-years WMAP data release is used to 
identify and analyse such regions. Two areas are selected with signal-to-noise ratio 
S/N < 2 and S/N < 3, covering ~ 16% and ~ 26% fraction of the sky, respectively. 
The polarization power spectra of these two areas are dominated by the sky signal 
on large angular scales (multipoles £ < 15), while the noise prevails on degree scales. 
Angular extrapolations show that the synchrotron emission competes with the CMBP 
i?-mode signal for tensor-to-scalar perturbation power ratio T/S = lO^'^-lO^^ at 
70-GHz in the 16% lowest emission sky {S/N < 2 area). These values worsen by a 
factor 5 in the S/N < 3 region. The novelty is that our estimates regard the whole 
lowest emission regions and outline a contamination better than that of the whole 
high Galactic latitude sky found by the WMAP team {T/S > 0.3). Such regions 
allow T/S ^ 10~^ to be measured directly which approximately corresponds to the 
limit imposed by using a sky coverage of 15%. This opens interesting perspectives to 
investigate the inflationary model space in lowest emission regions. 

Key words: cosmology: cosmic microwave background - polarization - radio con- 
tinuum: ISM - cosmology: diffuse radiation - radiation mechanisms; non-thermal. 



1 INTRODUCTION 

The Cosmic Microwave Background Polarization (CMBP) 
allows the study of the first stages of the Universe and 
is one of the hot topics of cosmology. One of its compo- 
nents, the _B-mode, is sensitive to the primordial Gravi- 
tational Wave (GW) background left by inflation giving a 
way to in vestigate the physics o f the very early Universe 
(e.g. iKam ionkowski & Kosowskvl ll998l) . In fact, its power 
spectrum on degree scales has a linear dependence on the 
tensor-to-scalar perturbation power ratio T/S, which mea- 
sures the amount o f primordial GW (e.g. iBovle et alJl2006t 
iKinnev et ai]l2006ll . In combination with parameters mea- 
sured by the CMB Temperature spectrum (the scalar pertur- 
bation spectral index n and its running dn/dlnk), it allows 
us to distinguish among the numerous models of inflation. 

The B-mode signal is fain t and only u pper limits 
of T/S have been set so far. ISeliak et al.l il200ell and 
iMartin fc Ring eval ( 2006), using different analyses, find con- 
sistent results of T/S < 0.22 and T/S < 0.21 (95% C.L.), 
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respectively. On the other hand, the plethora of inflation 
models does not help constrain T/S. This can vary by or- 
ders of magnitude and ca n be even smaller than 10~*, al- 
though js^l^^i] i2006l) show that only models with high 
fine tuning degree can feature T/S < 10~^ . Correspond- 
ingly, the peak B-mode signal can vary from ~ 200 nK of 
the present upper limits down to values even smaller than 
3 nK for T/S < 10^*. 

The weakness of the CMB B-mode makes it easily con- 
taminated by foreground emissions from both the Galaxy 
and extragalactic sources. The study of astrophysical fore- 
grounds is thus crucial to set the capability of CMBP exper- 
iments to investigate the early Universe and help set which 
part of the infiationary model space is accessible. 

At frequencies lower than 60-70 GHz, the most rele- 
vant contaminant is expected to be the Galactic synchrotron 
emission. Recent results based on 3-years Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) data at 22.8-GHz show 
that even the high Galactic latit udes are normally strongly 
contaminated (iPaee et alj|200^ . These authors use about 
75 percent of the sky (all high latitudes but strong large 
local structures like the Northern Galactic Spur) and find 
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Figure 1. Top. P olarized intensity map at 22.8-GHz of the 
WMAP experiment iPaefe et al.|2(^)y ^ .Data a re averaged on ~ 4° 
scale using the HEALPix lG6rsk^taI .|2005^ pixelation with pa- 
rameter Nside = 16 (resolution r4). The map is in Galactic co- 
ordinates with the Galactic centre in the middle. The coordinate 
grid is 15° spaced both in longitude and latitude. Lowest emission 
regions are the blue spots at mid-high Galactic latitudes out of 
the large SNR loops. Mid and bottom. Lowest emission regions in 
the 22.8-GHz WMAP polarization map: the 4°-pixels of the two 
regions with S/N < 2 (mid) and S/N < 3 (bottom) are blanked 
(white) . 



that at 60-70 GHz the synchrotron emission competes with 
the cosmic B-mode signal even for models with T/S — 0.3- 
0.5, which are already disfavoured by the present upper lim- 
its. Similar results are obtained by iLa Porta et al.l (l200(il) 
through the analysis of the 1.4 GHz sky at latitude > 20° 
using the DRAG survey data JWolleben et al.ll2006l) . 



iLa Porta et alJ i200d) also analyse areas^ with smaller 
extension (for a total of about 10% of the sky) and lower 
emission, finding that the Galactic signal at 70-GHz is com- 
parable to the GMB B-mode in case of T/S ~ 0.1. Even 
though lower, this signal would largely contaminate the 
GMB B-mode for most of the T/S values presently allowed. 
However, these areas are not in the lowest emission regions 
visible in the 22.8-GHz WMAP polarization map. 

Information about the lowest emission regions are 
available, instead, from three independent small areas ob- 
served at 1.4 and 2.3-GHz jGarretti et alj l2005bl l2006l : 
iBernardTe t al. 2006|). The synchrotron emission here actu- 
ally looks weaker, competing with the GMB B-mode only 
for cosmological mod els with T/S = 10"^-10"^ at 70-GHz 
ijGarretti et al.ll200(ih . 

Such values would give more chances to investigate a 
large part of the inflationary model space, opening interest- 
ing perspective for the study of the inflation physics. How- 
ever, the areas observed so far cover just few square degrees 
(~ 20-deg^ in total) and represent few sparse samples of 
the lowest emission regions. A real knowledge of the typical 
Galactic synchrotron contribution of such regions is still an 
open question, and requires wider areas to be analysed. 

This paper reports the first study carried out to iden- 
tify and analyse the areas with the lowest synchrotron emis- 
sion at microwave wavelengths, which allow negligible Fara- 
day rotation effects ( still appreciable at ~ 1-GHz; e.g. see 
iGarretti etHI l2005a^ . We use the 22.8-GHz polarization 
map of the 3-years WMAP data release to select them (Sec- 
tion I^J, and compute their polarized angular power spec- 
tra (Section Finally, we estimate the contamination of 
the GMBP by Galactic synchrotron in the GMBP frequency 
window (70-90 GHz) and discuss implications for the de- 
tection of the cosmic B-mode in lowest emission regions 
(Section^. 



2 LOWEST EMISSION REGIONS SELECTION 

The 3-years WMAP data release^ has made available the 
first polarization all-sky map at m icrowave frequencies, 
namely 22.8-GHz JPage et al.ll200d) . The sensitivity does 
not allow to detect the signal out of the Galactic plane at 
the nominal resolution (FWHM ~ 1°). The high Galactic 
latitudes have low signal-to-noise ratios (S/N), except in 
large local structures like the big Supernova Remnant (SNR) 
loops. The situation improves once data are averaged on 
pixels of ~ 4° (top panel of Fig ure as presented by the 
WMAP team jPaee et alJl2006l) . At this resolution the sig- 
nal appears even at high Galactic latitudes and allows us to 
identify the lowest emission regions in the sky (blue areas 
in top panel of Figure 0, which, however, still have signal 
competing with the noise. 

We quantitatively select the lowest emission regions us- 
ing the 4° polarized intensity map. After accounting for the 
noise bias dUvaniker fc Landeckeil [20021) . pixels with S/N 
smaller than a given threshold have been selected. Namely, 

1 Three areas defined as follows: A) 180° < I < 276°, 6 > 45°; 
B) 193° < I < 228°, b < -45°; C) 65° < Z < 180°, b > 45°; I and 
b are Galactic longitude and latitude. 
^ http://lambda.gsfc.nasa.gov/ 
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we have considered two values: S/N < 2 and S/N < 3. The 
two resulting regions are shown in Figure 

These regions are not contiguous and located at mid- 
high Galactic latitudes, mainly in the third Galactic quad- 
rant, with significant extensions in the second one. It is 
worth noting that the Galactic caps are not included, having 
too high emission. 

The sky fractions covered are 16.2% {S/N < 2) and 
25.7% {S/N < 3), respectively. Although smaller than half 
a sky, these two regions have significant extensions and could 
still effectively be used for CMBP B-mode studies. 



3 POWER SPECTRUM ANALYSIS 

We compute the E- and _B-mode power spectra of the po- 
larized emission in the two selected areas. These are the 
quantities predicted by cosmological models and allow a di- 
rect comparison with the CMBP signal. 

To account for the irregular sky coverage we use the 
method based on two-point correlation function s of the 
Stokes parameters Q and U described by Sb arra et alJ 
i2003h . Correlation functions are estimated on the Q and 
U maps of the selected regions as 

c''{e) = xaj x = Q,u (i) 

where Xi is the pixel i content of map X, i and j identify 
pixel pairs at distance 6. Data are binned with pixel-size 
resolution. Power spectra Cf'^ are obtained by integration 

Cf+Cf-"^Wt / [C'^{e)Fi,e2{e)+C"{e)F2,e2{S)]sm{e)de{2) 
Jo 

C^+Cf-"^We / [C"{e)Fi^e2{e)+C'^{e)F2,i2{S)]sm{0)dO{3) 
Jo 

where the f unctio ns Fi^im and F2,em are described by 
IZaldarriaeal (Il998l) . Cf and Cf are the noise spectra, 
and We is the pixel window function accounting for pixel 
smearing effects. The WMAP map smoothed on 1° pixels 
(IfEALPix parameter Nside — 64) has been used to com- 
pute correlation functions and spectra, instead of the 4° one 
used for the selection. This allows us to reduce the effects of 
the pixel window function and provides reliable computation 
up to ^ ~ 100 of the B-mode peak. 

Figure |21 shows B-mode spectra in the two selected 
areas. These are steep on large angular scales, with the 
S/N < 2 case having less power than S/N < 3, as expected 
according to selection criteria. A flattening instead occurs 
on small angular scales. This behaviour indicates the signal 
prevails on large scales {£ < 10 and ^ < 15 for S/N — 2 and 
3, respectively), while the noise dominates on smallest ones. 
To check if the flat component is compatible with noise, we 
performed Montecarlo simulations. We have generated 100 
noise map realizations assuming the sensitivity map pro- 
vided in the WMAP data package. Then we compute their 
spectra using only the pixels of the two selected regions. For 
both the two cases, we find that the average spectrum is flat 
with normalization = 7.5 x 10^^ fj.K^ and 7.3 x 10"^ fiK^ 
for S/N < 2 and S/N < 3, respectively, which are com- 
patible with the level of our spectra at large £ (Figure I^J. 
In addition, the spectra of the two S/N cases converge to 
the same value in the flat range, which again suggests the 
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Figure 2. B-mode spectra of the WMAP 22.8-GHz polariza- 
tion map in the two regions with S/N < 2 (solid) and S/N < 3 
(dashed). To avoid confusion, error bars are reported for the 
S/N < 2 case only. The mean noise spectra of the Montecarlo 
simulations are also reported (dot— dashed and triple-dot-dashed 
for S/N < 2 and S/N < 3, respectively). 



Table 1. Mean value of (^£{1 + 1)/ {2n)Cf') - see text. Fit pa- 
rameters of S-mode spectra for S/N < 3 are also reported. The 
f-range used for the estimates is indicated. 



S/N 


£-range Afg [^K^] 






2 
3 


[4, 12] 

[4, 24] 0.30 ± 0.06 


-2.4 ±0.3 


2.8 ±0.3 
7.62 ±0.14 



spectra are noise dominated there. We use the mean noise 
spectra attained with the Montecarlo simulations to account 
for the noise bias in all the next evaluations. 
We tried to fit a power law 




to the spectra. Only the case S/N < 3 allows meaningful 
results (Table 0, because only few points (three) can be ef- 
fectively used for S/N < 2, after the noise bias is subtracted 
for. To provide an estimate of the mean emission also for the 
S/N < 2 case, we compute the mean value of the quantity 
£{£ + l)/(27r)Cf over the usable ^-range, which is £ ^ 12. 
The result is reported in Table Q along with the same es- 
timate for the S/N < 3 area. It is worth noting that there 
is a factor ~ 3 in spectrum (~ \/3 in signal) between the 
mean emissions of the two analysed regions. Moreover, the 
emission we found in the S/N < 2 region is a factor ~ 10 
weaker than t hat measured in t he ~ 75% sky fraction in the 
same £-range iPa.ee et alJl2006l) . 



4 IMPLICATION FOR CMBP B-MODE AND 
DISCUSSION 

The contamination of the CMB B-mode signal can be esti- 
mated by extrapolating the spectra of Sect. El to 70- GHz, a 
frequency in the range where the combined con tribution of 
synchrotron and dust is believed to be minimum dPaee et alJ 
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I2OO6I: ICarretti et'ai] I2OO6II . It is worth noting that typi- 
cal values of rotation measure at high Galactic latitudes 
{RM ~ 10-20 rad m~^) generate polarization angle rota- 
tions of 6-12 arcmin at 22.8-GHz. Thus, Faraday rotation 
effects can be considered negligible, allowing safe frequency 
extrapolations. 

Our spectra are noise dominated at the multipole I the 
CMB B-mode peaks at [l ~ 90) and a direct extrapolation 
of them would provide just upper limits. Then, we adopt an- 
other approach: we use the spectra measured in the ^-range 
where they are signal-dominated and perform angular ex- 
trapolations to the scale of interest. The mean emission is 
that provided by the last column of Table Q As for the ex- 
trapolation, the analysis of the 1.4-GHz DRAO survey data 
shows that the angular behaviour of Galactic synchrotron 
spectra is well represented by a power law 

Ci oc (5) 

with slope 13 varying in the range [-2.5, -3. 01 for I = 2-300 
iLa Porta et al. 2006; Buri gana et aljl2006l) . This is consis- 
tent with the value of /3 = -2.6 quoted by t he WMAP 
team for the whole 75% of the sky they consider jPage et alJ 
[200^. We use as slope both the two edge values of the 
j3 range (-2.5 and -3.0) and start extrapolations from the 
middle of the ^-range used to evaluate the mean spectrum, 
namely io = 8, 14 for S/N < 2, < 3, respectively. 

Frequency extrapolations are performed assuming the 
brightness Temperature of the Galactic sync hrotron follows 
a pow er law T^y'^"'^ oc i/", with a = -3.1 iBernardi et alJ 

Results are reported in Figure |21 For the S/N < 2 
case, which concerns ~ 16% of the sky, the contamination 
range covered by the two angular extrapolations at ^ ~ 90 
competes with the CMBP _B-mode for models with T/S = 
[3 X 10"'^, 1 X 10~^]. The comparison with previous mea- 
surements shows these values are significantly lower than 
the mean contamination estimated by the WMAP team for 
most of high Galactic latitudes. On the other hand, our re- 
sults are similar to the contamination estimat ed in the small 
sky p atches observed in low emission regions dCarretti et alJ 
I2OO6I and references therein) giving a valuable consistency of 
results in regions selected with similar criteria (areas in low- 
est emission regions) . Although with the uncertainties of an 
angular scale extrapolation, our results extend to the whole 
lowest emission regions the promising indications obtained 
in those few small sample areas covering just few square de- 
grees, and strongly support the possibilities to access small 
T/S values in the lowest emission part of the sky. 

The use of frequency slope a = — 3.1 can be conside red 
as a conservative approach. In fact. lHinshaw et al.l i200(fl fit 
the behaviour of the synchrotron total intensity, finding out 
a steepening of the slope at the highest frequencies of the 
WMAP range. In particular, they fit a = —3.5 in the range 
23-60 GHz. With such a slope the synchrotron contamina- 
tion at 70 GHz further reduces by a factor 2.5 with respect 
to a = -3.1 and even T/S ~ [1 x 10"^, 3 x 10"^] becomes 
accessible (yellow shaded area of Figure |HJ . 

The situation is slightly worse for S/N < 3, as expected, 
since higher emission pixels are included in the analysis. 
However, T/S = [2 x 10~^ 4 x 10"^] looks still accessi- 
ble even in this larger region, assuming a = —3.1. These 
are values better than that for all the high Galactic lati- 
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Figure 3. B— mode power spectra of the Galactic synchrotron 
emission estimated at 70-GHz in the regions with S/N < 2 
(green shaded area) and S/N < 3 (gray shaded area). The 
two angular extrapolations with /3 = —2.5 and —3.0 are re- 
ported (upper and lower dashed line of shaded areas). A slope 
a = —3.1 is assumed for the frequency extrapolation. The case 
of a slope a = —3.5 is also reported for the S/N < 2 region 
(yellow shaded). For comparison, the plot shows CMB spectra 
for three different values of r = T/S, the synchrotron emission 
detected in the two low e mission areas in the target fields of the 
BOOMERanG CCg^ p. . . ICarretti et al. 2005b) and BaR-SPOrt 
experiments (C^^^^, Icarret^^^al![p006 | ) . and the up per limit 
found in the DASI fields (Cg.^p iBernardi et alj|2006l) . Finally, 
the general contamination at high Galactic latitude, as estimated 
by the WMAP team using 74.3% of t he sky is also repo rted: the 
spectrum at 22.8-GHz is considered iPage et al.ll200a) and ex- 
trapolated to 70-GHz. 



tudes (WM AP results on ^ 75% o f the sky) , or in the areas 
analysed bv lLa Porta et al.l (|200(f)- 

For the lowest emission areas, the minimum combined 
synchrotron-dust contamination has been found to be some- 
where between 70 and 90- GHz I Carretti et'ai] 120061') . It is 
thus worthy estimating the synchrotron contamination of 
our best case {S/N < 2) even at 90-GHz. A further factor 
~ 4 is gained with respect to 70-GHz (in thermodynamic 
temperature), and the Galactic signal competes with mod- 
els with T/S ~ [1 X 10"^ 3 X 10~^]. 

Despite the significant general contamination of high 
Galactic latitudes, the results obtained by our analysis de- 
pict a situation with better and interesting conditions in the 
lowest emission regions. The possibility to reach T/S values 
as low as 10~"^-10~^ looks now extended from few small 
good sky patches to an area representing ~ 16% of the sky, 
that could be the right place for deep CMBP observations 
looking at the B-mode. 

To limit observations in small regions, however, im- 
poses intrinsic constraints on the minimum detectable T/S, 
mainly because of the leakage from E- into the weaker 
_B-mode. In fact, lAmarie et al.l l|2"005) find that an all-sky 
survey would allow a detection of T/S with a theoretical 
sensitivity limit of A{T/S) = 1.5 x 10~^, that becomes 
A{T/S) = 3.2 X 10"^ when 70% of the sky is available, and 
A{T/S) = 10"^ for 15% (3-a C.L.). It is worth noting that 
the class of infiationary models with minimal fine-tuning 
have T/S values ranging within 10~^ and 10~^ iBovle et alJ 
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l2006h . for which a 15% sky portion would be large enough 
for the first detection of the tensor CMBP component. 

In spite of these limitation, the minimum T/S value 
detectable in 15% of the sky almost matches the limit im- 
posed by the foregrounds we have found in lowest emission 
regions. Although the uncertainties because of the angular 
extrapolations we need to apply, the search for the _B-mode 
in a sky portion of such a size could thus be a good trade-off 
between intrinsic and foregrounds limits. 

In addition, the weakness of the B-mode signal makes 
already a challenge to dete ct the signal for T/S = 0.1 w ith 
the present technology (e.g. lCortiglioni fc Carrettill2006l) . It 
is likely that significant technological improvements will be 
necessary before cosmologists can face an all-sky mapping 
mission with a sensitivity able to match the intrinsic limit 
of A{T/S) ~ 10~^. Therefore, an experiment aiming at de- 
tecting the _B-mode in a smaller region (10-15% of the sky) 
can be a valuable intermediate step that would allow us to 
probe inflation models down to T/S = 10~^. 
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